Roto-vibrational spectrum and Wigner crystallization in two-electron parabolic 

quantum dots 
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We provide a quantitative determination of the crystallization onset for two electrons in a parabolic 
two-dimensional confinement. This system is shown to be well described by a roto-vibrational model, 
Wigner crystallization occurring when the rotational motion gets decoupled from the vibrational 
one. The Wigner molecule thus formed is characterized by its moment of inertia and by the corre- 
sponding sequence of rotational excited states. The role of a vertical magnetic field is also considered. 
Additional support to the analysis is given by the Hartree-Fock phase diagram for the ground state 
and by the random-phase approximation for the moment of inertia and vibron excitations. 

PACS numbers: 73.21. La, 73.21.-b 



I. INTRODUCTION 

For a low enough electron density, Wigneri predicted 
that electrons should localize creating an ordered spa- 
tial structure, thenceforth named a Wigner crystal, that 
breaks the complete translational symmetry of the ho- 
mogeneous electron gas (also see Ref. 0). Indeed, the 
formation of the Wigner crystal was observed in two- 
dimensional (2D) distributions of electrons on the sur- 
face of liquid helium.^ A phase transition, induced by 
the magnetic field, from an electron liquid to a crystalline 
structure has also been reported for a 2D electron plasma 
at a GaAs/AlGaAs heterojunctioni4 

The existence of different phases in quantum dots, 
where a few electrons are confined into a small space, 
has become a topical subject in mesoscopic physics (see, 
for a recent review, Ref. In fact, the high controlla- 
bility of quantum dots suggests that these systems could 
provide an attractive opportunity to achieve the appro- 
priate conditions for localized states. It is precisely to 
stress this controllability that the names artificial atoms 
and quantum dots have been coined. 

There is a general persuasion that the Wigner crys- 
tallization in quantum dots, whose localized states are 
referred to as Wigner molecules, should occur at signif- 
icantly larger densities than in the 2D bulk. It is based 
on the argument that in quantum dots potential-energy 
contributions can easily exceed the kinetic terms and, 
therefore, electronic motion can be effectively quenched 
by manipulating the external confinement and/or an ap- 
plied magnetic field. As for the homogeneous gas, one 
would expect that in crystallized states the kinetic energy 
is solely that of the vibrational zero-point motion of the 
electrons around their mean positions, much smaller than 
the interaction (potential) energy. Various approaches 
including ab initio calculations within diffusion and path 
integral Monte Carlo methods, Hartree-Fock and spin- 
density functional methods etc have been applied to an- 
alyze the onset of the crystallizationn^ However, a non- 



ambiguous theoretical result that would justify the above 
conjecture for a zero magnetic field is lacking. The case 
with an intense magnetic field is better understood since 
the magnetic field induces an edge reconstruction, be- 
ginning with the appearance of localized vortices on the 
outer region, that ultimately propagates to all the dot 
for very high JS'sj^ 

In the simpler case of a two-electron 2D quantum 
dot at zero magnetic field, Yannouleas and Landman^ 
pointed out that the excited-state energies of this system 
closely follow the rotor sequence when the repulsion-to- 
confinement ratio, as given by the Wigner parameter Rw, 
is large enough (^-^200). This was shown to be a proof of 
the crystallization of the two electrons on fixed positions 
in a reference frame which is rotating. Quite remark- 
ably, the hypothesized rotating Wigner molecule fulfills 
at the same time the strict symmetry conditions of quan- 
tum mechanics -circularity in this case- and the obvious 
preference for opposite positions when repulsion is large 
enough. This is a major difference from the above men- 
tioned bulk case where a Hamiltonian symmetry (trans- 
lation) is broken by the crystallized state. For Wigner 
molecules, symmetries are preserved in the laboratory 
frame and one must consider an intrinsic (rotating) frame 
to 'see' the underlying deformation. A similar situation 
is found for particular states of two-electron atoms that 
have been much investigated in physical chemistry (we 
address the reader to the review paper by BerrjsSi). For 
the two-electron quantum dot, however, the crystalliza- 
tion condition from Ref. 0, Rw ~ 200, looks disappoint- 
ing since it seems unrealistic to achieve such a value ex- 
perimentally. 

Although the exact ground-state wave function of the 
two-electron artificial atom can be obtained, at least nu- 
merically, it may seem paradoxical that one also needs 
the excited states in order to ascertain the existence of 
a crystallization. In fact, this inability to disentangle 
the system's intrinsic structure from its full wave func- 
tion in a clear way can be taken as a weakness of the 
ab initio, symmetry preserving, approaches. In general. 
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even in cases when the exact ground- and excited-state 
wave functions and energies are known, an intrinsic de- 
formation can only be inferred by comparing with the 
resuh of simpler models in which either symmetries are 
relaxed or the intrinsic structure is imposed. A clear ex- 
ample of the former approach is given by the unrestricted 
Hartree-Fock (HF) method for the ground state^-^" fol- 
lowed by the random-phase approximation (RPA) for 
excitations;^^ Conversely, the roto-vibrational model of 
Wendler et alm^ for two electrons in a ring could be in- 
cluded in the latter category. 

One should be aware that when symmetries are re- 
laxed, as in the Hartree-Fock approach, artifacts or un- 
physical properties may appear. In a recent contribution 
Reusch and Grabertii discussed the validity of the lat- 
ter, drawing special attention to the caution with which 
one must take Hartree-Fock predictions on symmetry 
breaking, in agreement with the results presented below. 
Therefore, a complete physical understanding requires 
both exact results and model solutions. This way the 
system's intrinsic deformations are physically understood 
while, at the same time, artifacts can be safely discarded. 
A paradigmatic case where the proposed analysis can 
be performed is given by the two-electron 2D parabolic 
quantum dot. The separation of center-of-mass and rel- 
ative coordinates along with the circular symmetry re- 
striction allows the reduction of the multidimensional 
Schrodinger equation to just a radial one, easily solvable 
numerically. On the other hand, the Hartree-Fock and 
RPA solutions without any symmetry restriction can also 
be obtained. A most convenient basis for this latter cal- 
culations is given by the Fock-Darwin orbitals in terms 
of which one can analytically develop much of the re- 
quired algebra. It is our aim in this work to determine the 
crystallization onset of two-electron parabolic dots by re- 
course to the three different approaches referred to above; 
namely, (a) an analytical roto-vibrational model, (b) a 
numerical solution of the Schrodinger equation, and (c) 
symmetry unrestricted Hartree-Fock and random-phase 
approximations. 

Hereafter, we refer to the solution of the Schrodinger 
equation for the two-electron parabolic dot as the exact 
solution. It should be pointed out that, as shown by 
Taut,^** this Schrodinger equation is analytically solvable 
only for particular confinement /interaction strengths. 
For general values of this quantity a numerical treatment 
is required. As mentioned above, the most straightfor- 
ward one is an integration of the radial equationii^ but, 
nevertheless, other methods such as diagonalization in 
a basisiSiii and the Monte Carlo method-'^^-^^ have also 
been applied. One of us has used the so-called oscilla- 
tor representation method, perturbatively treating the 
residual interaction, to derive analytical expressions for 
the energy levels 

The paper is organized as follows. Section II introduces 
the magneto-parabolic units that allow one to trace the 
evolution of ground and excited states of artificial atoms 
at various conditions. An analytical roto-vibrational 



model for the two-electron parabolic quantum dot is de- 
scribed in Sec. HI. Section IV provides details of our nu- 
merical calculation of exact solutions and compares these 
solutions with those of the roto-vibrational model. Sec- 
tion V analyzes the reliability of the Hartree-Fock and 
RPA results for the present system. A short summary is 
finally drawn in Sec. VI. 



II. MAGNETO-PARABOLIC UNIT SYSTEM 

We consider two electrons with a Coulomb interaction 
between them. The electrons move in the xy plane where 
a circular parabolic confinement induces the formation 
of an electron island. The system is also subject to an 
external magnetic field applied in the vertical direction 
{z). The full Hamiltonian thus reads 
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In Eq. (Q, m, K, and g* are the electron's effective mass, 
the dielectric constant, and the effective gyromagnetic 
factor, respectively, and we have used planar polar coor- 



dinates (r 



■y )■ The two contributions within the 



square brackets are, respectively, the generalized kinetic 
energy in terms of the vector potential A, and the exter- 
nal confinement. Within the so-called symmetric gauge 
one has A{x,y) = B /2{—y,x). The next contribution is 
the Coulomb repulsion and, finally, the last term is the 
Zeeman energy involving the total spin operator ^2 and 
the universal Bohr's magneton = e7i/2TOeC. 

It is well known that in the chosen gauge the mag- 
netic field contributions can be recast into the form 
of an effective parabolic confinement of frequency SI — 
■^/o^Q + i^c/^j where — eB/mc is the cyclotron fre- 
quency, and an additional angular- momentum-dependent 
term {hujc/'2,)iz (cf. ri,efs. l2ll|22 ). The magneto-parabolic 
units (mpu) we shall use consist of taking hft as the en- 
ergy unit and £fi = ^Jh/{mQ) as the length unit. In 
addition, one also imposes fi as angular momentum unit 
which, in turn, fixes the time unit Tfi — l/Sl. Summa- 
rizing in the standard abuse of notation we may write 
fi — Q = io, = 1 mpu. This is a natural choice for 
magneto-parabolic confinements and it allows one to ex- 
press the spatial part of the Hamiltonian in terms of only 
two adimensional parameters, namely. 
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Note that i?,„p and Wmp give, respectively, the ratios 
of Coulomb interaction strength and cyclotron frequency 
to effective confinement. In the absence of a magnetic 
field Rmp coincides with the so-called Wigner parameter 
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FIG. 1: Equivalence between the mpu pair of adimensional 
paramaters (K-mp,Wmp) and the physical values of the ex- 
ternal parabolic strength and magnetic field {hu!o,B). The 
bulk GaAs effective mas m — 0.067me and dielectric constant 
K — 12.4 have been assumed. 



Wrap = 2 that corresponds to a zero confinement wq = 0. 
We also mention that Reusch and Grabert used these 
parameters in their recent Hartree-Fock calculations 
Within the mpu system the Hamiltonian reads 
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where m* is the adimensional ratio of effective to bare 
electron mass, i.e., m* — m/rrie. The passage from the 
mpu system, with a given Rmp and Wmp, to physical 
units requires the knowledge of the effective mass m and 
the dielectric constant k. More specifically, with fixed m 
and K one can invert Eqs. 10) and (PJ for the effective 
confinement and cyclotron loc frequencies or, equiva- 
Icntly, for the external confinement ojq and the magnetic 
field B. In terms of these results the physical values of 
the energy {Ml) and length units are readily found. 

Figure ^ shows the equivalence between the adimen- 
sional parameters {Rmp , Wmp ) of the mpu Hamiltonian 
and the physical values of {Tiluq, B) for the case of a GaAs 
host semiconductor having m = 0.067me and n = 12 A. 
The advantage of working with the mpu system becomes 
obvious when realizing that, with unit redefinition, the 
same numerical results can be equally applied to a variety 
of confinements, magnetic fields, and material parame- 
ters (effective mass and dielectric constants). Therefore, 
the model acquires a certain degree of universality. We 
also expect that quantum dot properties such as scaling 



laws or phase diagrams will be better displayed in terms 
of the adimensional mpu parameters. 



III. A ROTO- VIBRATIONAL MODEL 

Taking Eq. ^ as a starting point and introducing the 
standard center of mass {R, 8) and relative (r, 6) coor- 
dinates it is well known that the Hamiltonian separates 
and, therefore, that the wave function factorizes. The 
center of mass (CM) problem is that of a single particle 
in a harmonic potential and magnetic field, having an 
analytical solution in terms of Fock-Darwin orbitals and 
energies e^]^' =2N+\M\ + 1 + MWmp/2, with (iV, M) 
the radial and angular momentum center-of-mass quan- 
tum numbers. Focusing next on the relative problem, 
one introduces the wave function e™^w„m(r)/-yf hav- 
ing good £z angular momentum (m), and an additional 
quantum number n whose meaning will be clarified be- 
low. The equation for the unknown u„m('') reads 
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■ rnWmp/'i in terms 
of the relative- motion energy Snm' and the mpu param- 



where we have defined £nm "* 
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Rw of Ref. [3- Also note that Wmp has a maximal value ^^^^ 



Equation ^ will be the basis of our roto-vibrational 
model. Note that it resembles a Schrodinger one- 
dimcnsional equation with an effective potential 

(6) 
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that includes the rotational motion term ~ rr? jr^ char- 
acterized by the angular momentum quantum number 
TO. We can expect a rigid-rotor behaviour if Veff{r) has 
a deep minimum at a particular value r = tq. When this 
occurs the situation resembles that of diatomic molecules 
like H2, where the potential well for nuclear motion is 
described by the Morse potential (see Ref. "^s) . In the 
present case the effective potential indeed has a minimum 
although it is in general rather shallow. This property is 
responsible for the coupling between rotation and vibra- 
tion or, equivalently, for the floppiness of the rotating 
molecule mentioned in Ref. 0- 

The minimum condition on Veff {r) yields the rotor ra- 
dius from 

ro_ _ Rmp _ 2(to^ 
2 
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Neglecting the third contribution on the left-hand-side, 
an assumption that will be valid for large enough tq, one 
finds the asymptotic law ro ~ {2Rmpy^^- Now, expand- 
ing to second order around tq we approximate 

l/4\ „ ,2 
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a result that, when substituted into Eq. JSJ for the 
round parentheses, leads immediately to the analytical 
predictiooSi 

?(rm) _ 1 2 , R^np - I / i 

The result embodied by Eq. (jS)) has a clear physi- 
cal interpretation. It contains a rotor-like contribution, 
~ 771^/ {2J'), with a moment of inertia given hy J' — 
mpu, and a vibrational one characterized by a quan- 
tum number n. The vibrational frequency LUyn, = \fkf\i 
(// = 1/2) is given by the last square-root factor. Similar 
to atomic molecules, there is roto- vibrational coupling, 
since the vibration frequency depends on m and, in ad- 
dition, centrifugal distortion since Tq also depends on to. 
For large enough values of i?™^, implying large rg and 
therefore small average densities, the centrifugal distor- 
tion disappears and one has tq ~ i^Rm-pY^"^ f^^' ^ vrCs. 
In this limit rotational terms become negligible, as well as 
roto-vibrational ones. Thus, Eq. reduces to a simple 
TO-independent asymptotic expression 

sl™^ = ^^'™^l + ^^(" + • (10) 

When adding the magnetic field, the roto-vibrational 
energy becomes 

^ \ 2) ' ' 

in agreement with the expectations from Ref. jl^i for two 
interacting electrons in a quantum ring^ It is worth 
stressing that since the Wmp dependence only amounts to 
an energy shift of the Eq. lO eigenvalue, the radial func- 
tion Unmiy) does not depend on Wrap- Therefore, one 
may conclude that the roto-vibrational properties are not 
affected by magnetic fields, for a fixed (n, to) state. Of 
course, since the energy shift varies for different states, 
the magnetic field will modify the ordering of energy lev- 
els. For instance, the level crossings as a function of 
Wrap will cause the ground state to have a nonvanishing 
TO-value. This actually explains the buildup of increasing 
permanent currents in the dot's ground state. 

The results from this section will be validated by com- 
paring with the exact ones below. The roto-vibrational 
model presented here allows one to determine the crystal- 
lization onset from the criterion that rotation and vibra- 
tion motions decouple when intrinsic-frame electron lo- 
calization sets in. Conversely, when the coupling is strong 
the system could be represented by either a vibrating ro- 
tor or a rotating vibrator and, therefore, the situation can 
not be clearly resolved. It is also worth stressing that the 
roto-vibrational model describes all possible excitations 




FIG. 2: (Color online) Upper panel: Ground state energy in 
mpu's for the two electron dot. Lower panel: Ground state 
relative angular momentum. Note that since M = 0, relative 
and total angular momentum coincide. Even-m (odd-m) re- 
gions correspond to singlet (triplet) ground states. Results for 
1.9 < Wrap < 2 are not shown due to the excessively large vari- 
ations of the computed quantities in this region. The dotted 
lines separate in each domain with a given m the crystallized 
(above) from the non crystallized (below) phases using the 
criterion of roto-vibrational coupling below 3 % (See subsect. 
IV.B). 



of the relative-motion problem. For this particular sys- 
tem, this amounts to a description of all the excitations 
since the center-of-mass and spin degrees of freedom can 
be analytically integrated out. 
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r (mpu) 



FIG. 3: (Color online) Upper panel: Efltective potential de- 
fined by Eq. @ for m = and two diflFerent values of Rmp- 
The relative-motion eigenenergy ioo is given in each case 
by the horizontal line. Lower panel: radial eigenfunction 
uoo{'r)/^/r for the same two potentials of the upper panel. 



from the exponential decay from some (large) r onwards. 
In spite of this possible difference in asymptotic behavior, 
the two methods (matching and nodecounting) provide 
to a high accuracy the same eigenvalue. 

With the above numerical methods an exploration of a 
part of the Rmp — Wmp plane has been performed. Figure 
|21summarizes the results for ground-state energy and an- 
gular momentum. Note that the ground state always has 
(iV, M) = (0, 0) for the CM quantum numbers and that, 
because of symmetry, even (odd) m states are associated 
with singlet (triplet) total spin. Clear singlet-triplet os- 
cillations, as studied in the literature (cf . Ref. llTj) , are 
seen in the lower panel of Fig. [21 It is also worth men- 
tioning that the ground-state energy contour lines are 
piece- wise linear as a result of the simple dependence on 
the Wmp parameter. Indeed, the energy is totally Wmp- 
independent for m = and a fixed Rmp- 

Figure 131 displays the radial wave functions uoo(f)/\A' 
for two different values of Rmp as well as the correspond- 
ing effective potentials Veg defined by Eq. We note 
that, in agreement with the discussion of the preceding 
section, when increasing Rmp the effective-potential min- 
imum moves outwards and it effectively binds the lower 
states to its neighborhood. When this occurs the radial 
probability is strongly quenched at small r (lower panel) 
and the scenario indeed resembles the familiar one from 
the physics of diatomic molecules. A more detailed com- 
parison of the roto-vibrational model with the exact re- 
sults will be done in the next subsection. 



B. Crystallization criterion 



IV. ONSET OF WIGNER CRYSTALLIZATION 
A. Exact solutions 

We have solved the radial equation for the relative 
problem, Eq. lO, numerically. Several standard ap- 
proaches to eigenvalue problems with boundary condi- 
tions can be used for this purpose. Specifically, we have 
applied the so-called matching method where one in- 
tegrates with the Numerov or Runge-Kutta algorithms 
from the origin r = outwards assuming an 
behaviour. Additionally, imposing an exponential decay 
law exp (— r^/4) for large r, inwards integration is per- 
formed and the two solutions are required to match at an 
intermediate r-point. To ascertain the numerical result 
for the eigenvalue, as a control we also used the method of 
'node-counting', where only outwards integration is per- 
formed and the eigenvalue is found from the condition 
that Unm {t) increases by one the number of radial nodes 
from the required value when the energy exceeds the cor- 
rect eigenvalue by an infinitesimal amount. In principle, 
the node-counting method assures the correct boundary 
condition at r — > oo automatically. However, in practice, 
there always remains a small difference between the nu- 
merical and exact eigenvalues, responsible for a deviation 



The validity of the roto-vibrational model presented in 
Sec. Ill is proved by the results of Fig. 01 As shown in 
the upper panel, the error of the analytical prediction for 
Enm is important only when Rmp is small. As a matter of 
fact, for Rmp > 2 the discrepancy for ego is always below 
2%, even with the asymptotic expression Eq. H1U|I . Note 
that the relative errors slightly increase with increasing 
n (dashed curves), although the analytical approxima- 
tion can still be considered quite good. The analogous 
comparison for to > (not shown) yields smaller relative 
errors than those for m = 0. The lower panel of Fig. 
^analyzes the excitation energies as measured from eooj 
again showing an excellent agreement between the exact 
results and the ones obtained with the analytical model, 
with small deviations only at small Rmp- Based on the 

asymptotic expression, Eq. (|10|) . we have chosen 1/Rmp 
as an independent variable in order to better display the 
linear behavior associated with the rigid rotor at large 

Rmp - 

We shall rely on the high accuracy of the roto- 
vibrational model to provide a quantitative measure of 
the crystallization onset. Our criterion will be the fol- 
lowing: the two-electron parabolic system is assumed to 
be crystallized as a rotating Wigner molecule when the 
roto-vibrational coupling falls below a given percentage 
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FIG. 4: (Color online) Upper panel: inm exact energies and 
the analytical prediction of the roto-vibrational model (thin 
solid curves). The dashed lines give the relative error in per- 
centage (right scale) with bottom to top curves corresponding 
to n = 0, 1, ... 4, respectively. Lower panel: Excitation ener- 
gies of the (0, m) states, i.e, eom — £oo from the exact (data) 
and analytical model (curves). 



(typically chosen as 2 or 3 %). The roto-vibrational cou- 
pling is defined as 



7(m) = 100 



Wvibini; Rr, 



V3 



V3 



(12) 



where ujmb is the square-root factor in Eq. and -v/S 



is the limit of this quantity for i?„ 



CG. Using a 3% 



condition the crystallization onset for each angular mo- 
mentum is given by the Rmp value where each m curve 
of Fig. Renters the shaded region. Note that the crystal- 
lization onset moves towards higher values as the angu- 
lar momentum is increased, reflecting the property that 
roto-vibrational coupling is stronger for high-m states. 
It is also worth mentioning that, since not all m-states 
are simultaneously crystallized, in practice the rotational 
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FIG. 5: (Color online) Roto-vibrational coupling as defined in 
Eq. The shaded region indicates the crystallized phase 

with the criterion 7 < 3. 



hands will gradually degrade for increasing m's as a con- 
sequence of the roto-vibrational coupling. In agreement 
with the bulk gas situation, the crystallized states with 
the proposed criterion are characterized by having a po- 
tential energy that largely exceeds the kinetic one, as can 
be easily checked from Eq. I^. 

The crystallization properties of m and — m states are 
identical since one can easily check that Unmii", Rmp) = 
Un-m{r; Rmp)- Therefore, taking into account the varia- 
tions in ground state angular momentum we can draw the 



rap 



mp 



boundaries for Wigner crystallization in the R 
plane, i.e., the crystallization phase diagram (see the 
lower panel of Fig. Of course, if instead of a 3% 
threshold for roto-vibrational decoupling one chooses a 
different value the crystallization onset will vary, al- 
though as shown in Fig.|Sl for low m's the crystallization 
is not crucially dependent on the precise percentage in 
the range 2 to 4%. Actually, it should be more appropri- 
ate to speak of crystallization onset for a given percentage 
of roto-vibrational decoupling than of an absolute value. 

It should be noted that the onset of Wigner crystalliza- 
tion may be studied by means of the conditional proba- 
bility distribution (CPD, cf. Rcfs. ll0ll2l for finding one 
electron at positon ri given that the second electron is at 
position Y2 ■ As soon as the interaction is switched on, the 
CPD, which for a system with only two electrons is just 
the modulus squared of the wave-function |\E'(ri, r2)P, 
exhibits the formation of a molecular-like state (see Fig. 
EJ. It is difficult, however, to use this measure alone 
as a conclusive evidence of the formation of the Wigner 
molecule in two-electron dots. The difficulty arises from 
the fact that even weak interactions, for which we should 
not expect a crystallization, yield the formation of a hole 
around the electron at Y2 (the correlation hole) and a 
maximum at ri = — r2. Indeed, as seen from Fig.|^ the 
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FIG. 6: (Color online) Contour plots of the ground state two- 
electron wave function j'I'(ri,r2)P for the fixed value of r2 
shown by a solid symbol. A value of Wra-p = and the 
displayed Rm-p have been used for the different panels. For 
Rm-p = 0, r2 has been arbitrarily fixed on the x axis, while 
for Rmp > it has been placed at the distance inferred from 
the asympotic law x-i = ro/2 (see text). 

results for Rmp ~ 1 hint at a molecular state even though 
according to our analysis the roto- vibrational coupling is 
still strong and this state is not associated with a crys- 
talhzcd phase (see Fig.lSJ. When R^^p increases the de- 
pletion of the CPD around the fixed electron becomes 
much stronger, which is in qualitative agreement with 
the crystallization trend from our analysis. 

V. HARTREE-FOCK AND RPA APPROACHES 

In this section we discuss the results obtained within 
the symmetry unrestricted Hartree-Fock method for the 
ground state and the corresponding RPA for excitations. 
We recall that the HP and RPA approaches were devised 
for the analysis of many-body systems. Therefore, their 
application to a two-electron quantum dot is merely an 
exploratory approach to the qualitative features of the 
Wigner crystallization rather than a quantitative descrip- 
tion of the above exact results. 

We have solved the HP problem in the Pock-Darwin 
basis which diagonalizes the square bracket in Eq. 
namely {\arj)] a = l,...M;ri =t,i}, where a labels 
the orbital part and rj the spin. Our basis has been 
optimized such a way that we consider the 70 lowest 
Pock-Darwin states, of the non-interacting energy level 
scheme, for a chosen value of the magnetic field. An 
arbitrary single-particle orbital |i) is then expanded as 

\i) = J2arj^ari In the chosen basis, the HP equa- 

tions are written as a system of nonlinear eigenvalue 
equations for the matrix of B coefficients (see, for in- 
stance, details in Ref . ITl|) . 

We have imposed good Sz HP orbitals leaving totally 
unspecified the remaining spatial symmetries. Note that 
the Slater determinant built with these single-particle or- 
bitals will be an eigenstate of the total Sz operator but 
not, in general, of S^. This, as we shall see in the results, 
is intimately connected with the prediction of broken cir- 
cular symmetry. 

In Pig. [71 we show the HP phase diagram in the 



Rmp — Wmp plane (lower panel) and the corresponding 
total energies (upper panel) . The total energies resemble 
those obtained in the exact treatment, with approximate 
piecewise linear regions between orbital angular momen- 
tum transition lines. As expected the actual values at a 
given point in the diagram lie slightly above the corre- 
sponding exact results. In the lower panel different gray 
regions reflect a measur^^ of the deviation from circular- 
ity of the ground state density, with the lightest intensity 
corresponding to a circular (nonbroken symmetry) solu- 
tion and more intense gray levels to noncircular (broken 
symmetry) results. The contour lines show the total or- 
bital angular momentum L^. Regions I and III are of 
circular symmetry and for them Lz has a good quantum 
number, taking the values and -1, respectively. Por 
the rest of the diagram (regions II, IV and V) the con- 
tour lines only indicate the expectation value of Lz but, 
since the mean field in not circularly symmetric, this is 
no longer a good quantum number. The dotted curve 
separates the states having total spin projection 5*^=0 
(below) from Sz = I (above). 

Por Sz ^ configurations, the HP method predicts 
a broken-symmetry solution when Rmp exceeds a value 
~ 1, somehow below the onset of crystallization obtained 
in Sec. IV {Rmp — 2.2). In this region (II) the HP solution 
is indeed a mixture of singlet and triplet states, as can be 
verified by computing the total spin dispersion AS^. The 
corresponding spatial density is built from two opposite 
and localized single electron orbitals. It is instructive to 
compare the HP solutions in region II with those obtained 
using a total-spin conserving ansatz 

•i'sngiri,r2,rii,ri2) = (t){ri)(j){r2) Xsngim,V2) (13) 

-xxtrpivum) (14) 

for singlet (5'sng) and triplet (^'trp) states, where 
A[(f>i{ri),<j)2{j^2)] denotes the antisymmetrized product 
of the two orbitals (j>i and (f>2 while the x's are the well 
known singlet and triplet spin states. 

As a sample result the total energies obtained for 
{Wmp, Rmp) = (0,2) are E,,act = 3.720, Ehf = 4.034, 
Esng = 4.185 and Etrp = 4.168; i.e., by requiring total 
conservation the mean field energy raises considerably. In 
addition, while the HP solution breaks circular symme- 
try, both ^E'sng and "i'trp yield circular densities because of 
the spatial dependence of the 0's in Eqs. H13|l and H14|l . 
We have also checked that these results are equivalent 
to those obtained using the Lipkin-Nogami projection 
method^^'^^ for the effective Hamiltonian Hcs — 7Y — AS^ 
in order to restore symmetry approximately. The 
above ansatz for states with good total spin are examples 
of the use of constraints in mean field approaches, which 
necessarily raise the energy above the mean field mini- 
mum. One could also project the symmetry- unrestricted 
HP orbitals as discussed in Ref. |33- In the latter case, 
however, the wave function is no longer a single Slater 
determinant but rather a sum of few determinants of 



8 




J ■ ^ ^ ■ . ■ ^ ^ 1 

1 2 



W 

mp 

FIG. 7: (Color online) Upper panel: HF ground state energy 
in mpu's for the two electron dot. Lower panel: Regions (I 
to V) of the HF ground state phase diagram. The dotted 
curve separates total Sz = (below) and Sz = 1 (above) 
configurations. The gray scale denotes the space symmetry 
of the HF solution, from circular (light gray) to strongly 
deformed (dark gray) configurations. Labeled contour lines 
display (Lz). Note that in the broken-symmetry regions (H, 
IV and V) this latter quantity varies continuously between 
the integer boundaries. Results for Wmp > 1.8 are not shown 
due to the difficulty in determining the single-particle basis 
when the Fock-Darwin states become quasi-degenerate. 



the corresponding symmetry operator. Therefore, the 
ground state energy with the restored symmetry is no 
longer bound by the mean-field minimum and thus can 
be closer to the absolute minimum imposed by the vari- 
ational principle. 

With the above results, we conclude that the lowest HF 
solution in region II requires a simultaneous breaking of 
the spin and space symmetries. Taking into account the 
results of the preceding sections we can say that singlet- 
triplet mixing in region II is an artifact of the HF solu- 
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FIG. 8: (Color online) Moment of inertia computed in the 
RPA (circles) at Wm^ = as a function of the adimensional 
parameter -Rmp. Solid lines show the evolution of the corre- 
sponding values in the analytical model of Sec. HI for different 
m-states. Each line starts at the crytallization onset for the 
corresponding angular momentum, according to the criterion 
of Sec. IV. B. The dashed line represents the asymptotic value 
J — r-o/2 taking ro « (2i?mp)^''^, to which all solid lines 
converge at very high values of Rmy 



tion. However, space symmetry breaking in this region is 
a true physical mechanism indicating an intrinsic struc- 
tural change in the exact wave function, as supported 
by the roto-vibrational model discussed above. We be- 
lieve this peculiar combination of artifact and physics is 
due to the smallness of the configuration space for a two- 
electron system. We must also point out that at high 
Rraf and/or Wmf (regions IV and V) the HF prediction 
fails to match the results of the analysis given in Sec. IV. 

In contrast to the exact results of Fig. |21 where a re- 
gion (although narrow) with m = —2 corresponding to a 
singlet state appears at large Wmp, the HF approxima- 
tion predicts only one {Sz = 0) ^ {Sz = 1) transition 
as Wmp increases. This can be understood as an overes- 
timation of the exchange energy in the HF model which 
tends to favor spin alignment whenever orbital overlap- 
ping occurs, as it does in region HI with circular orbitals, 
as well as in regions IV and V with two-lobed orbitals. 

We consider next the results obtained by solving the 
RPA for excitations. As discussed in Ref. (llj the RPA 
determines the moment of inertia associated with the col- 
lective rotation of a deformed HF structure (see Eq. (18) 
in Ref. 11). Figure|Hlpresents the evolution with Rmp, at 
Wmp = 0, of the RPA moment of inertia Jrpa (circles). 
For comparison, the values computed through the solu- 
tion of Eq. {T)) of the roto-vibrational model, J — rg/2, 
are also shown (solid lines). Each J'-line starts at the 
crystallization onset for the corresponding angular mo- 
mentum using the criterion of Sec. IV. B. Note that Jr,pa 
remains null until the HF solution breaks rotational sym- 
metry at Rmp — 1, from there on it reasonably agrees 
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with the exact values, somehow averaging the exact re- 
sults for different m's. The molecule stretching, yield- 
ing larger tq's (J^'s) as m increases, is obviously outside 
the RPA. All J' values slowly converge to a common re- 
sult with increasing Rmp, i-C, to an exact rigid-rotor be- 
havior. For comparison, the dashed line represents the 
asymptotic value corresponding to tq « (2i?mj,)-'^/^. The 
nice qualitative agreement between Jrpa and J is giv- 
ing additional support to the above conclusion that space 
symmetry breaking in region II of the HF phase diagram 
indicates a genuine physical effect and, thus, it also sup- 
ports the overall picture of a rotating Wigner molecule. 

Although the RPA restores circularity on the deformed 
HF mean fieldfii associated with the operator, it is 
not able to restore the symmetry related to S^, since 
the latter one is a two-body operator which is beyond 
the RPA treatment of broken symmetries. A side effect 
of the HF spin artifact in region II is that when spin- 
flip bosonic pairs are included in the RPA, generalizing 
our previous calculation^ the correlation energy in this 
symmetry-broken phase is badly overestimated, i.e., it is 
between two and three times larger than the exact value. 
This does not occur, however, in the circularly symmetric 
regions. 

Focusing now on the RPA vibron states, we must also 
distinguish excitations associated with spin and space 
(charge) degrees of freedom. While the RPA qualita- 
tively describes all excitations in regions I and III, it fails 
for spin excitations in phase II. Obviously, this is due to 
the HF spin artifact in this region. It also fails for space 
excitations of phases IV and V. Finally, we end this sec- 
tion by pointing out that, in spite of the spin artifact 
of phase II, in pure triplet phases RPA reproduces the 
exact spin precession frequencies known from the the- 
ory of magnetic resonancei^ That is, a pure spin-flip 
state (precessional mode) is expected at the Larmor en- 
ergy Tiuoc = 9* IJ-bB = \g*m*Wmp mpu. Indeed, within 
the RPA, i.e., in the quasi-boson approximation ) ^^i'^^ one 
finds 

[n,ot]^nLocO+ , (15) 

with the vibron operator for the Larmor mode 

Sx + iSy 



(16) 



In Eq. I|16() (Sz) is the HF expectation value of the Sz op- 
erator. In fact, the Larmor mode at hujc appears when- 
ever the ground state has (Sz) ^ and it is normally the 
lowest excitation of the system. 



dots subject to an external magnetic field. The analysis 
has been done in terms of magneto-parabolic units and 
the associated parameters {Rmp, Wmp) that give, respec- 
tively, the ratios of Coulomb interaction strength and cy- 
clotron frequency to effective confinement. The ground 
state calculations are summarized in a phase diagram 
that can be equally applied to a variety of confinements, 
magnetic fields and material parameters. 

We suggested an analytical model for the interpre- 
tation of the exact results, including roto-vibrational 
coupling and centrifugal distortion (molecule stretch- 
ing). Within this roto-vibrational model we proposed a 
criterion to determine the onset of Wigner crystalliza- 
tion based on the decoupling of rotational and vibra- 
tional motions. For a 3%-decoupling threshold we found 
that Wigner crystallization appears, for zero-angular- 
momentum states, when Rmp exceeds a value ~ 2. States 
with larger m's crystallize at higher Rmp values. In agree- 
ment with the homogeneous gas situation the potential 
energy of the crystallized states is much larger than the 
kinetic energy, the latter one being solely due to the vi- 
brational zero-point motion of the electrons. 

The HF calculations predict that crystallization for 
Sz = occurs when Rmp > 1, the new phase (II) being 
in an artificial mixture of singlet and triplet spin states. 
The space symmetry breaking in phase II is a genuine 
physical effect but the spin mixture is an artifact due to 
the smallness of the configuration space for a two-electron 
system. Other HF symmetry-broken phases (IV and V) 
do not agree with the exact results. The RPA moment 
of inertia qualitatively agrees with the result from the 
roto-vibrational model, although the molecule centrifu- 
gal distortion is missed. On the other hand, the RPA 
produces reliable results for space (charge) excitations 
in regions I, II and III, as well as for spin excitations 
when the HF solutions possess good and Sz quantum 
numbers (regions I and III). We would expect a broader 
applicability of the many-body theories (HF-I-RPA) for 
larger systems. Work along this line is in progress. In 
conclusion, the combined use of exact and model calcu- 
lations allowed us to ascertain the existence of a rotating 
Wigner molecule in a two-electron dot for relatively large 
electron densities or, equivalently, small i?,„j, parameters. 
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